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SUMMARY 


A theoretical model is developed to predict losses in radial 
inflow turbine nozzles, T-he analysis is presented in two parts. 

The first one evaluates the losses which occur across the vaned 
region of the nozzle, while the second part deals with the losses 
which take place in the vaneless field. 

In the vaned region, equations are derived to relate the 
losses to the boundary layer characteristic parameters on the 
vane surfaces and over the two end walls. The effects of vanes 
geometry, flow conditions and boundary layer characteristic 
parameters, at exit from the nozzle channel, on the level of 
losses are investigated. Results indicate that the portion of 
the losses incurred due to the end walls boundary layer may be 
significant especially when these boundary layers are characterized 
by a strong cross flow component. 

In the vaneless field, the governing equations for the flow 
are formulated using the assumptions of the first order boundary 
layer theory. The viscous losses are evaluated through the 
introduction of a wall shear stress that takes into account the 
effects of the three dimensional end walls boundary layer. The 
resulting governing equations are solved numerically and the 
effects of some operating conditions and the nozzle geometry are 
studied. 

It is concluded that the losses in the vane less region may 
be looked at as a secondary factor in the determination of the 
overall efficiency of the turbine nozzle assembly. 
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INTRODUCTION 


In recent studies dealing with the flow patterns in radial 
inflow turbines, the flow is assumed to be isentropic with 
a viscous loss coefficient introduced to account for actual 
flow effects. The studies of Ref. [1, 2] were based on the 
assumption that the viscous losses within the flow passages 
are proportional to the average kinetic energy in the passage. 

A proportionality constant was chosen such that the viscous 
losses matched the values determined from the experimental 
results. Other studies [3], estimate losses using approxi- 
mations based upon flat plate or pipe-flow analogies. Although 
in the past, approximate and semi-empirical solutions for the 
losses have been adequate for engineering purposes, the present 
sophistication of fluid machinery demands a closer view of 
the problem. 

The real flow in the radial rype machines is very complex 
because it is viscous, unsteady and three dimensional. An 
insight into the real flow processes can be achieved by breaking 
the complex flow pattern into several regions. Existing flow 
theories can then be used to obtain a solution to the flow 
together with some reasonable simplifying assumptions in each 
of these regions. Following this approach, and referring to 
Figure 1, the flow field in a radial turbine nozzle can be 
divided into the following subdomains: 

a. The passage between the blades, in which the flow 
is partially guided. 

b. The throat region where the flow is directed from 
the closed channel conditions to the vane outlet 
free conditions . 

c. The region just at the exit of the channel, which 
will be referred to as the zone of rapid adjustment. 

In this region, a considerable portion of the 
overall stator losses take place. 

d. The vaneless space region. 

In the present work a theoretical model is developed to 
predict the losses in the flow regions mentioned above. The 
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investigation is presented in two parts. The first part deals 
with the flow losses in the vaned regions and the second part 
with the losses in vaneless field. In the first part/ the 
friction losses which result from the boundary layers develop- 
ment along the vane surfaces and over the end walls are 
considered. Mixing losses, due to the nonuniform flow 
conditions at exit from the vaned region, are also taken 
into account. The analysis establishes relations that provide 
the description of the flow properties at the entrance to the 
vaneless field, in terms of the boundary layer characteristic 
parameters over the flow surfaces. From these relations 
the various factors influencing the losses in the vaned region 
are determined, and their relative effects are evaluated. 

In the second part of this study, the viscous flow in 
the vaneless field is analyzed. The influence of the three 
dimensional end wall boundary layers on the steady nozzle flow 
is taken into consideration. The effect of the nozzle 
geometry and the flow characteristics on the losses in the 
vaneless region are determined. 

Under the stipulation of the developed model," the overall 
losses between the inlet of the radial nozzles and the impeller 
periphery may be considered as the sum of friction and mixing 
losses in the vaned region, plus the friction losses in the 
vaneless field. For an actual expansion process, the mechanism 
of all these losses are coupled. The model adopted however, 
provides a simple guide for better understanding and estimating 
radial nozzles loss data. 
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PART I 


FRICTION AND MIXING LOSSES IN THE VANED REGION 

ANALYSIS 

The purpose of the radial turbine nozzles is to convert the 
pressure energy into flow kinetics energy. In the vaned region, 
a torgue is imposed on the flow by means of the vanes following 
the scroll and hence, the angular momentum of the fluid changes 
more rapidly than in the vaneless field. The design objective 
remains the same, however, to achieve the highest possible 
efficiency. The thin untwisted vanes forming the nozzle channels, 
are arranged around the circumferential direction, with equal 
spacing and orientation (Fig. 1) . In order to maintain a 
sufficiently high vane chord Reynolds number, the length to 
width ratio of the passages between the vanes ranges between 
1 and 3, and the vanes aspect ratio is often less than 1. 

As the flow passes through the vaned region, it suffers a 
loss in total pressure and kinetic energy. This may be attributed 
to the decrease in the flow momentum within the boundary layers, 
which are developed over the flow surfaces, and to their mixing 
with the main flow. It can be seen from Figure 2, which is taken 
from Reference [4] , that the magnitude and the direction of the 
flow velocity, as well as the other flow properties are nonuniform 
in both the axial and the tangential directions at nozzle inlet. 
These nonuniform flow conditions which constitute an additional 
source of the losses are caused by the three dimensional flow 
behavior in the scroll. 

Another possible mechanism that contributes to the losses 
is the secondary flow, which tends to move the boundary layer 
material across the channel from the pressure side of a vane to 
the suction side of the adjacent vane, as shown in Figure 3. 

To obtain an exact solution of the real flow in the vaned 
region, with all the aforementioned effects included, would be 
impractical. These factors can be examined separately to 
determine how they affect the losses in the radial turbines. 

The work presented here is a study of the losses, which result 
from boundary layer development along uhe vanes surfaces and 
over the end walls. 
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Referring tc Figure 1, the total pressure of the flow entering 
the nozzle at radius r , is P t . As it leaves the nozzle at r-^, 
there is a variation in the flow properties in both the circumfer- 
ential and the axial directions due to the axial directions due 
to the boundary layer development over the end walls and the vane 
surfaces. It is assumed that due to the mixing effects, the 
properties of the air stream are homogeneous at station 2, down- 
stream from the vanes exit. The loss coefficients, expressed 
in terms of the total pressures at stations 0 and 2, will hence 
represent the overall losses across the vaned region, resulting 
from the friction and mixing mechanisms. The pressure and enthalpy 
loss coefficients are given respectively by the following 
relations : 


Y 


1 >P t2 //P to 

1 “ P?/ P to 




(P to/ P t2 

* P to^ p 2* 


) (y-D/y 
(y-D/y 


1 


1 


( 1 ) 


In order to determine these coefficients in terms of the main 
flow properties and the boundary layer parameters at station 1, 
the equations governing the flow motion between stations 1 and 2 
are derived in the following sections. 

The Governing Equations 

Consider a control volume with an axial depth, 2b, equal to 
the nozzle depth, that extends circumferentially between two 
consecutive mid-channel stream surfaces as shown in Figure 4. 

The equations of conservation of mass, energy, angular momentum 
and linear momentum in the Y-Y direction are written under the 
following assumptions: 

a. The total temperature is constant during the expansion 
and the mixing processes. 

b. The static pressure is independent of z and 0 at 
station 1. 

c. The flow pattern is similar in all the flow channels, 
between two neighboring vanes of angular spacing 
equal to 2ir/Z. 

The Conservation of Mass 

For steady flow conditions, the following continuity equa- 
tion applies for the control volume of Figure 4: 
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(2) 


+b +tt/Z +b +ir/Z 

f f (pVrcosa)-, d0dz = / / p 9 V„r_cosa 9 dedz 

-b -ir/Z X -b **tr/Z 

If the left hand side of the above equation is expressed in terms 
of the boundary layer characteristic parameters at station 1 
by using equation (A24) of Appendix A, the following relation 
is obtained: 

( p V) m r-^cosa^tl-6 -6 te -X8 ] = P2 V 2 r 2 GOSa 2 , (3) 

In the above equation, the subscript refers to the conditions 
outside the boundary layer regions at exit from the nozzle 
channels, a-, is the angle between the velocity vector and the 
inward radial direction, 8 is the total nondimensional dis- 
placement thickness over the vane surfaces at station 1 and 
6 denotes the trailing edge blocking factor. 8 is the 
streamwise total nondimensional momentum thickness over the two 
end walls at station 1, and X is a parameter given by Equation 
(A23) that depends on the end wall boundary layer profiles. 

Conservation of Angular Momentum 

If the equation of conservation of angular momentum is 
written for the same control volume, the following expression 
is obtained: 

+b +ir/Z 2 2 + ' !r / z 2 2 

/ / p 2 r 2 V 2 cosa 2 sina 2 dedz / S (pr V cosasina) -^dedz = 0 

-b -ir/Z -b -ir/Z ^ 

Using Equation (A23) of Appendix A, Equation (3) can be written 
as follows: 

2 22 . 2 * *** 
P 2 V 2 GO Sa 2 sina 2 r 2 ~(pV ) ra cosa-^sina^ r-^El-5 ] - 0 (5) 

* 

In the above equation, 6 is the total nondimensional momentum 
thickness over the vane surfaces and n is a parameter to be 
determined from the end wall boundary layer velocity profile 

at station 1, using Eq. (A31) . 
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Conservation of Linear Momentum 


Referring to Fig. 4 the momentum equation in the Y--Y 
direction is expressed as follows : 


d8dz 


+b 

+ 1T/Z 
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f 

/ 

[pV cos 

-b 

-tt/Z 



+b 

+ TT/Z 


+ / 
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- / 
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-b 

-ir/Z 


+b +ir/Z 
/ S 

-b -ir/Z 


( 6 ) 


The first and second terms in the right hand side of Equation (6) 

represent the component of pressure forces acting on the 

« 

surfaces CD and FE of Fig. 4, in the Y-Y direction. The next 
term Q, is the resultant of the pressure forces on the two 
mid channel stream surface CE and DF in the same direction. 

The contribution of the body forces is represented by C and the 
flow turning angle is equal to S . When the first term on the left 
hand side of Equation (6) is expressed in terms of boundary 
layer characteristic parameters according to Equation (A41) , 
Equation (6) reduces to: 


2p 


2tt 


2 V 2 r 2 cosa 2 cos (a 2 -B) sin — - (pV ) m cos r^ 

{ siwz (1 —re**, — (a* + s te + e*) cos 1} 


2pi r i 


sm 


- 2p r 


2*2 cos6 sin i + ft + it + h 


(7) 


where r is a parameter that depends on the end wall boundary 
layer profiles as given by Eq. (A42) . 
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The Energy Equation 


Assuming no heat or work exchange between the flow in the 
control volume and the surroundings and using the equation of 
state for a perfect gas we obtain: 




(p t J -. 



( 8 ) 


Where the subscript, t, refers to the total conditions. Assuming 
the viscosity effects to be negligible in the free stream 
between stations 0 and 1, P^ will be equal to and the 

energy equation can be written as follows: 


<P 

Tp 


tl J 


t2 



(9) 


Equations (3) , (5) , (7) and (9) will be used to calculate 
flow properties at station 2 in terms of the free stream and the 
boundary layer characteristic parameters at station 1. In order 
to carry out a solution the radius, r 2 , at which complete mixing 
has occurred is to be known. Moreover, the detailed flow 
conditions within the control volume are also needed in order 
to evaluate the quantities P, Q and the angle $ in Equation (7) . 


SUDDEN MIXING ANALYSIS 

The results obtained from the experimental investigation 
of Ref. [4] indicate that the nozzle wakes almost disappear 
near the zone of rapid adjustment, which is shown in Fig. 1. 
Assuming complete mixing to occur at a radius r 2 which is very 
close to r^ is therefore an authentic model of reality. 

For the limiting case of r 2 approaching r, , the flow 
deflection angle, the contributions of pressure forces acting 
on the mid-channel stream surfaces CE; DF as well as the body 
forces in the Y-Y direction will diminish. Consequently, the 
terms Q, C and the angle 3 in Equation (7) will tend to zero, 
and the governing equations are simplified to the following 
form: 
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(10) 


P t2 COSct x Pl M 1 (1-5 -5 te “XB ') 

p = * * 

to cosa 2 P 2 M 2 


* * **, * . * * ** 
M 2 sma 2 [l-S -6 te ~xe ] = smo^l-5 -6 -e “ n6 ^ 


( 11 ) 


* *2 2 ** ^ + ^+- e + 8 ^ z " 

M~ cosa, {M, cos^a, [1 -ra — 1 , 

z. 1 1 , 7T 2y 


tan | 


i + e 1 ' 1 - ^ m ! 2 > > 

= M* oo S0Ii a- S *-6 te -XB**)lM; 2 cos 2 a 2 + |±i (1 - M* 2 ) 1 


(12) 


where 

* V 1 
M = C— -- 
1 l V 


cr. 


* n 

p i= 


* V 2 

jyi = _ 

2 V 


cr. 


(“ — ) 2 
p t 2 


(13) 


Equations (11) and (12) can be solved simultaneously to determine 
. . * 

the axi symmetric flow properties M 2 and « 2 , providing that the 
flow conditions and the boundary layer characteristics are 
defined at the nozzle exit. Once M 0 is obtained the density 
ratio p 2 could be calculated from the following identity: 


* 



(1 


T-1 

Y+l 


M.V" 1 


(14) 


If a similar expression for in terms of M 1 is used together 
with Equation (14) into Equation (10) , the ratio between the 
total pressures P t2 and is obtained. These pressure, ratios 
are substituted into Equation (1) to determine the overall loss 
coefficient Y, c, across the nozzle vanes. These loss coefficients 
‘“ill hence depend upon the nozzle geometry, the flow conditions 
M 17 a and the boundary layer characteristics 6,0,9 , r, 

1 and n at the nozzle channel exit. 
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RESULTS AND DISCUSSION 


The results of the present analysis show the effects of 

variation in nozzle geometry, flow conditions M, , a, , and the 

£ ^ *1* -L 

boundary layer characteristics 5,6,6 , r , X and tj on the 

level of losses across the vaned region. A series of figures 

in which the variation of the loss coefficients Y, ? with only 

two of the afore -mentioned variables are presented at specified 

values for the rest of the parameters. 

Figures 5, 6 and 7 represent the overall loss coefficient 

Y and ? at flow discharge angle, ct^, of 45°, 60° and 75°. Each 

figure demonstrates that the overall loss coefficients increase 

with increased nondimens ional vane surface boundary layer momentum 
* 

thickness, 8 , for different values of the streamwise nondimens ional 

* * 

momentum thickness of the end wall boundary layer, 6 . It can 

also be seen from these figures that the losses decrease with 

the decrease in the end wall momentum thickness 8 , with the 

* * 

minimum values corresponding to 8 = 0 . This lower limit 

represents the losses in a nozzle channel of infinite aspect 

* * * 

ratio. The nondimens ional momentum thicknesses 8 and 6 are 
defined by Equations (A22) and (A27) and are both evaluated at 
exit from the nozzle channels. 

✓ 

Equation (A27) shows that the nondimensional momentum 
* 

thickness, 6 , is inversely proportional to the cosine of the 
flow discharge angle a, . For a specified value of the vane 
surface boundary layer momentum thickness, 8 will increase 
with increased a^. Figures 5, 6 and 7 can therefore be used to 
show that the losses increase with increased flow discharge 
angles. This effect becomes more significant at higher values 
of a^. An accurate estimation of the exit flow angles is hence 
essential for loss evaluation. 

Figure 8 shows the effect of the streamwise shape parameter, 

H , on Y for the case of small cross: flows., It is: obvious; 

X 

the Y is a weak function of H , which is known not to exceed 2,. 

A 

for unseparated turbulent boundary layers [6] . For all values 
** 

of 8 presented, the loss coefficients were found to increase 
not more than 5% for a corresponding increase in H x from 1.2 to 2. 


REPRtJDl3CIBILITY OF THH 
OR PAGE IS POOR 
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Since the vane trailing edge thickness, t , is increased 

when cooling of the radial turbine is considered, it would be 

desirable to examine the effect of t on the overall losses 

e 

in some detail. This effect is investigated through the use 

of the parameter t -f-S +6 / ^J ' Z . The. result's are presented 

g s p r *i ^ 

in Figure 9 as the ratio between the loss coefficient, Y, at 

a finite value of t to Y, calculated for zero trailing edge 

e 

thickness. Figure 9 shows that the increase in the losses can 

be considerable for larger values of the trailing edge thickness. 

* 

Under normal operating conditions corresponding to 6 =0.02 

and 9 = 0.015 in a conventional turbine with trailing edge 

2ir Z 

thickness parameter, t +5 +5 /= — , of about 0.03 the 
c e S p 

trailing edge contribution to the losses is about 25 percent. 


2uZ 


If the parameter, t +6 +6 /■ 

0 S P JL I 


, is increased to 0.07, the 


portion of the losses due to the finite trailing edge is 
increased to 82 percent. 

The effect of compressibility on the overall loss coef- 
ficients are shown in Figures 10 and 11. It can be seen 
from Figure 10 that the enthalpy loss coefficient, Z , remains 
practically constant with increased free stream Mach number, 
for the different values of the boundary layer parameters 
considered. On the -other hand, the overall total pressure 
loss coefficient, Y, increases significantly with increased 
Mach number, Figure 11. Due to the large sensitivity of the 
pressure loss coefficient to Mach number variations, it is 
recommended to use the overall enghalpy loss coefficient, £, 
in describing the loss performance characteristics of the 
radial turbine nozzles. 

The effects of the end wall boundary layer cross flow on 

the losses are shown in Figures 12 and 13. Figure- 12 illustrates 

— * 

the variation of Y with 0 , for the different end wall 
velocity profiles A, B, C and D of Table A-l in Appendix A. 

Since it was found previously that the loss coefficients 
are weak functions of H it may be concluded, by examining 
Figure 12, that the cross flow profile' parameters K and L 
(see Equation a 9 and A10) have a significant effect on the 
losses. As an example, for a value of the nondimens ional 
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^ £ 

streamwise momentum thickness, 8 , of 0.03 a variation in the 

magnitude of losses- up to 25% can result by considering 
different profiles. 

Figure 13 shows the contribution of the cross flow to 
the losses which were calculated using the end wall velocity 
profile, B. It is clear from the figure that the minimum 
losses are obtained for the case of collateral end wall boundary 
layers e = 0. The losses increase appreciably for strong 
cross flow cases corresponding to higher values of e . 
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PART II 


FRICTION LOSSES ‘IN THE VANELESS FIELD 


The vaneless nozzle in a radial tu-rbine consists of a smooth 
walled axisymmetric passage with radial or conical surfaces. The 
passage depth may change with radius , but a little influence on 
both the pressure and velocity distribution is noticed, since the 
angular component of the fluid momentum is predominant. Basically 
the flow in the vaneless nozzle can be treated as a vortex motion, 
except for the skin friction which acts upon the flow boundaries. 

In overcoming the resistance along its path the flow loses part of 
its total pressure. The pressure loss occurs at the expense of 
both static and dynamic pressure [8] . 

Loss analysis for the turbine vaneless nozzle may be treated 
using the same methods developed for the compressor vaneless 
diffuser, provided that the area changes as well as the flow 
direction are taken into consideration. Conventional flow 
analyses for the vaneless diffuser [9, 10, 11] were based on the 
assumption of one dimensional flow. Total pressure losses were 
determined using a constant friction coefficient. In a recent 
study, Jansen [12] arrived at a more accurate evaluation for the 
friction coefficient and its variation along the flow path. 

His analysis however, requires much more detailed knowledge of 
the flow conditions than can actually be realized in practice. 

In the following study, the governing equations for the flow 
in the vaneless space are formulated using three dimensional end 
wall boundary layer uheory [6] . The resulting differential equations 
are solved numerically to obtain the flow properties within the 
vaneless field. The losses are evaluated taking into account both 
the wall friction and the momentum flux changes due to velocity- 
profile variation along the flow path. The effects of some 
parameters representing the operating conditions and nozzle 
geometry on the losses incurred within the vaneless field are 
also investigated. 
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THE THREE DIMENSIONAL FLOW IN THE END WALLS BOUNDARY LAYER 


Flow Pattern 

As the flow in the vane less field proceeds towards the 
impeller tip, the main streamlines turn as shown in Fig. 14. 

A balance is established in the main stream between the static 
pressure gradients and the centrifugal forces caused by the flow 
turning. The unbalance between the externally imposed pressure 
gradients and the centrifugal forces of low momentum end wall 
boundary layers, produces an acceleration component normal to 
the main stream direction in the boundary layer fluid. 
Consequently, the velocity in end wall boundary layers is 
different from the mainstream velocity, V, not only in magnitude 
but also in direction. The velocity profiles in this skewed 
boundary layer is shown in Fig. 15. Taylor [133 noticed that 
the polar plot of the velocity vector in a skewed boundary layer 
follows a triangular form such as shown in Fig. 15. Thus, 
any analysis of such boundary layer is usually carried out using 
different formulation in the inner and outer sublayers. 


The Inner Sublayer 

Referring to Fig. 15, the following linear relation between 
the cross flow velocity component, w, and the streamwise velocity 
component, v, holds good: 

w = v tan ill n 5 ^ 

y max '- LOJ 

The angle 4 , inax in the last equation is invariant within the 
collateral inner sublayer. Accordingly, it can be estimated 
at the wall using the following relation 


tan ij> 


max 


, , w 

Limit — 

z -*0 


e 


(16) 


Since the skin friction is in the same direction as the velocity 
vector at the wall, therefore it will deviate by the angle, 

^max' from the main flow direction. If the wall shear stress 
and the friction coefficient are denoted by and respectively 
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and if their components in the main flow direction by t and 
C fx' then: 

C f = C fx (1 + e2 ) V2 (17) 

where 


= 


1 2 
7 


and 


x 


'fx 


1 2 
7 P v 


(18) 


(19) 


Ludwieg and Tillman [14] deduced an accurate semi-empirical 
expression for the wall friction of two dimensional flows in 
the presence of favorable pressure gradients. Johnston [6] 
showed that in a skewed boundary' layer, the relation of 
Ref. [14] can be used to determine the component of the skin 
friction coefficient in the main flow direction according to 
the following relation: 

-0.268 

C fx = 0.246 [Exp t-1. 561 H^) ] R Qx (20) 


where H is the streamwise shape factor, and R. is the 
Reynolds number based on the streamwise momentum thickness . 

Referring to Fig. 16^ the components of shear stresses in 
the tangential and the radial directions are 


T e = 


t sin (a - ill _) 
o max 


(21) 


T = 

r 


t cos (ct - ill ) 
o r max 


( 22 ) 


When Equations (16), (17) and (18) are substituted in Equations 
(21) and (22) , the following relations are obtained: 


T 



pV 2 C 


fx 


[sina 


e cosa] 


(23) 
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( 24 ) 


1 2 . 

T r “ 2 1 pV C fx^ cosa + e s:Lna '* 

In the above equations, the values of a and e will depend 
on the main flow variables. The expressions which are used in 
evaluating these parameters are given in the outer sublayer 
analysis that follows. 

The Outer Sublayer 

The two velocity components in the non-collateral outer 
sublayer, namely the cross flow component, w, and the streamwise 
component, v, are related by the following equation, 

w = A (V - v) 


The parameter A depends on the main flow turning angle ($) , 
which is shown in Pig. 14, and is given by the following expression: 

A = - 2V 2 / % (25a) 

6 


where V is the main flow velocity. 

The tangent, e, of the limiting angle, ^ raax / can be expressed 
in terms of the parameter A and the streamwise friction coefficient, 
C^ x , according to the analysis of Ref. {6J as; 


£ = o.i i (1 


1/4 


-] - 1 


✓c 


fx 


(25b) 


Equations (20) , (25) and (25b) can be substituted into Equations 
(23) and (24) to predict the variation of shear components t q , x r 
along the flow path in terms of main flow variables. The details 
of the computation procedure used for this prediction is 
presented in the section dealing with the Numerical Solution of 
the Governing Equations . 
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CALCULATION OF THE MAIN FLOW PROPERTIES IN THE VANELESS FIELD 


Analysis 

The vaneless space geometry, together with the velocity 
tri-angles and the coordinate system are shown in Fig. 17. It is 
assumed that the flow enters the vaneless field at the radius, 
in a steady uniform axisymmetric pattern, with specified 
flow conditions. These inlet flow properties are determined 
from the vanes friction and mixing loss analysis presented in 
Part I of this study. 

As the flow passes through the vaneless field, there will 
be a variation in its properties in both radial and axial 
directions. The variation in the axial direction, z, results 
from the development of boundary layer over the two end walls. 
The static pressure distribution and the radial and tangential 
velocity components will be determined along the main flow path, 
using the assumptions of the first order boundary layer theory. 
An equivalent one dimensional flow will be considered in a 
nozzle with an effective depth of 2b e , which is equal to the 
actual nozzle depth, 2b, minus the displacement thicknesses 
of the two end wall boundary layers at any radius, r (see 
Fig. 18) . The components of the velocity, V, and the pressure, 
p, of the flow passing through this equivalent nozzle, vary 
only along the radial direction, and are independent of z and 
6 . The total pressure loss due to friction will be calculated 
using the values of and t q developed in the previous section 
and given by Equations (23) and (24) . 

Equations of Motion 

The control volume of Fig. 19 will be used in the following 
derivation. The magnitude of the velocity, V, is considered to 
increase as the flow proceeds inward in the vaneless nozzle. 
Furthermore, the flow will be assumed to be steady and adiabatic. 
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Conservation of Angular Momentum: 


Setting the net efflux of angular momentum through the 
control volume, equal to the external torque due to the shear 
force in the tangential direction, we obtain the' following 
relation : 

pb V cosa j— (r V sina) = x„r {26) 

e ar o 


Substituting for x. from Equation (23) into Equation 

9 


(24) , we get: 

V cosa V sina 


d C fx V 

+ V cosa (V sina) = r[sina - 


ECOSa] 

C 26 a) 


Conservation of Linear Momentum in Radial Direction : 

d pb V 2 sin 2 a 

-pb g V cosa ^ (V cosa) = b Q ^ (p) - x r - — C27) 

The left hand side of the above equation is equal to the net 
radial momentum flux. The first two terms on the right hand 
side are the radial components of the pressure forces and the 
wall friction forces, while the third term represents the 
centrifugal forces. 

Substituting for x r from Equation (24) into (27) and 
rearranging we obtain: 

-V cosa (V cosa) = — [cosa + esina]- ^ - S ^ ■ - - 

ar P ar r 

(27a) 


Conservation of Mass : 

From continuity considerations, the following differential 
equation is obtained: 
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I §£. + __1 

p dr V cosct 


0 


(28) 


d(V cosa) ,1,1 

35 + 5 + W JE~ 

e 


Equation of state : 

T-he pressure density temperature relation for a perfect 
gas is used in the following differential form: 


ldp = ldp ldT 

p dr p dr T dr 


(29) 


Energy Equation : 

Since the total temperature remains- constant during the 
expansion process in the vaneless field, the differential form 
of the energy equation may be written as : 


1 dT 
T dr 


Y-l 

2 

1 + 



2 dr 


M 


(30) 


Which, when combined with the definition of the Mach number 
reduces to: 


1_ dV^_ 

V 2 dr 


1 _ dM^ . 1 

M 2 dr 1 + ^ M 2 


(30a) 


Formulation of Relations Describing the Variation of Flow 
Properties Along the Flow Path ~~ - - - 

Equations (26) through (30) will be manipulated to obtain 
the following relations which are used to calculate the pressure, 
Mach Number and the main flow angle. 


Pressure Distribution 

Combining Equations (29) and (30a) , the following relation 
is obtained 

1 dp _ 1 dp 2 1 dM^ (21) 

p dr p dr j. + izi m 2 k 2 dr 

2 2 

Dividing Equation (27a) by V cos <s and using the identity. 



(32) 


EEPrtUCIBILrrY OF THE 
6KNOtfL pas® IB POOR 
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one gets an expression describing the variation of the tangential 
velocity component with radius as 


V cosa dr 


(V cosa) = 


1 

— 2 2 ~ 

yM cos a 


1 dp t "'fx (cosa + esina) 
p Hr “ ~2B~ ~ 


. 2 
tan a 


2 

cos a 


(33) 

Id 1 cl 

Substituting for — — from Equation (31) , y C os 'a dr ^ COSc< ) 

from Equation (33) into the continuity Equation (28), the 
variation of pressure with radius is deduced as : 


1 dp _ 
p dr 


M 


'fx (cosa + esinct) 


2 

sec a 


2 2 
M - sec 


2b 


2 

cos a 


- ( 


Y— 1 
~T~ 


i + 


M , 1 dM 2 
2) M 2 dr 


1_ 

b 


db 

€ 

dr 


(34) 


Mach Number Distribution 

In order to obtain the variation of Mach number along the 
flow path, an additional equation relating the pressure in terms 
of M and a is required. This relation is obtained by combining 
the equation of angular momentum (26a) with the equation of 
linear momentum (27a) . 


1 _ 1 dp = ^fx , JL v 1_ dV 2 

v 2 p dr 2b e cosa 2y 2 dr 


(35) 


Substituting Equations (30a) and (32) , into the last equation, 
it reduces to: 


1 dp _ yM 2 r , -1 1_ dM 2 C fx , 

p dr 2 x + Xli m 2 M 2 dr ^cosa 


(36) 


Mow, eliminating — between Equations (34) and (36) , the 

Mach number variation along the radius, r, is obtained from the 
resultant expression. 
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1 dM 
M dr 


1 

2 2 

M - sec a 


2 

• sec a 


2b 


fx 2 

seca (yM 


<2 -T 

tan a + etana)+ g~] 

e 

(37) 


Flow Direction 


The differential equation of flow direction may be developed 
using the following relation: 


tana = 


V sina 

V cos a 


(38) 


This last equation, can be written in the following differential 
form: 


1 d (tang) _ 1 . d (V sina) _ 1 d (V cosa) 

tana dr V sina dr V cosa 3r 


(38a) 


Substituting Eq. (26a) and (33) into' Equation (38a), the following 
relation is obtained: 


1 d (tana) 
tana dr 


sec a 1_ d]D _ sec a 

t m s p dr r 


(39) 


1 ( 3.0 

Substituting for — ^ from Equation (34) and rearranging, the 
equation describing the variation of flow direction with radius 
takes the following form: 


da = 
dr 


-tana 



2 


sec a 



C fx seca 
2b 

e 


[1 +(y-l)M + stana] + 


1 


db 
6 

dr 


(40) 


Nondimensional Form of the Governing Equations : 

Before solving the equations of motion, they will first be 
written in a nondimensional form- The nozzle dimensions will be 
normalized with respect to the inlet radius, and the pressures 

with respect to the total inlet pressure, p as follows: 
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Equations (34) , (37) and (40) can therefore be written in 
the following non dimensional form using the above expressions. 


1 dP 
P dR 


2 C 

yM r ^fx 

m 2 2 1 B 

yM - Sec a 


(cosct + esing) 
cos: 2 a 


2 

sec a 
r 


- ( 



1 + M 


1 dM 2 1 dB, 

zz s 3 r: j 

M 


(42) 


1 dM 
M dR 


1 + M 2 

_ —~2 
M - sec a 


2 

r sec a , rx 
+ — 


seca (yM 2 - tan 2 a + stana).+ g- ^jj|] 

(43) 


and 


da 

dR 


-tana 

2 2 
M - sec a 


.,2 Cj. seca - i j-d 

+ -^B (1 +( Tt- 1)M + * tanC * )+ 5 S 1 


(44) 


Furthermore , if the thicknesses of the boundary layers developed 
over the two end walls are taken to be equal, then, B, the depth 
ratio, and its derivative could be expressed as: 


2b - 26 


B = 


dB 

dR 


6 x 

"2d (=J) 
E2_ 

dR 


(45) 


where 


6^ is the streamwise displacement thickness. 
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PROCEDURE FOR THE NUMERICAL SOLUTION OF THE GOVERNING EQUATIONS 


In order to determine the flow properties within the vaneless 
field, the system of nonlinear differential Equations (43) , (.44) 

and the algebraic Equations (20), (2,5a) and (25b) will be solved. 

A fourth order Runge-Kutta algorithm is used to solve the 
simultaneous first order ^ differential Equations (43) and (44). 

The solution procedure consists of the following steps: 

1. The values of M, a, at the vaneless nozzle . inlet, are 
assumed to be known from the analysis given in Part I. 

2. At any radius R^, knowing 1YL, a^, C fx , e, B and ^ , 
the increments in M and a, over the radial distance (-AR) is 
computed using the following recursion formula: 

M i+1 - M. = - A| (Ki + 2K 2 + 2K 3 + K 4 > 

A P 

a i+l ” a i = ~ (KK 1 + 2KK 2 + 2KK 3 + ^ 4 ) 


(46) 

(47) 


where 


Kj = *<V M ± - =i> 

KK X = g(X it K i ) 

K 2 " £(E i - TT ' M i + 



Ctj 



kk 2 


g(R ± 


AR 
2 ' 


+ 


K 1 

T" 



K 3 = f < R i - ^ ' 

KK 3 = g(R. - A| 


M. 



M. 


r 


+ 


K 


2 


2 


a . 




(48) 
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K 4 = f(R ± - AR , M ± + K 3 


a ± + KK 3 ) 


KK 4 = g(R ± - A R, + K 3 , a ± + KK 3 ) 


Where f represents 


dM 

dR 


and g 


da 

3r 


and are evaluated according 


to Equations (.43) and (.44). respectively , using the known values of 

'clB 

C fx , A, e , B and ^ at the radius, R^. 


3. At the end of the interval, at the new -radial location 

dB 

R^ - AR, the quantities C^, A, z, b and ^ are evaluated using 
Equations (20), C25al and (258) respectively, as will be explained. 
The streamwise friction coefficient, , is calculated according to 
Equation (20) . The boundary layer characteristic parameters, R^ x and 
H , in this equation, can be evaluated using any method of 
boundary layer solution. In the present analysis, the method of 
Reference 15 was used for this purpose. Integrating Equation (25b) 
numerically, starting from the vaneles's; nozzle ih-let,. un-til any- 
desired radius, where the mean flow turning angle, g, is known, 
gives the corresponding value of A at R^ +1 . These values of C fx 
and A are then substituted into Equation (25c) , which is solved 
for e using an iterative procedure. 


4. Refined values for 


M i+1' “i+1 at (E i 


AR) are obtained 


by substituting the arithmetic average between R. and R. of the 

dB 1 

variables C fx , e, B and ^ in the recursion formula (46) and (47) . 


dB 


at R i+1 


5. The corresponding values for , A, e, B, 
are calculated as shown in Step 3, using the values of M. +1 , 
obtained from Step 4. 


6. Steps 4 and 5 are repeated until successive values of the 
computed flow variables are within the required accuracy. In the 
present study, an accuracy of 0.001 was maintained in the different 
variables, which was achieved after 2 iterations. 

7. The whole procedure (Steps 2 through 6) is repeated for 
small values of (-AR) starting from the nozzle inlet up to the 
impeller tip. The numerical results reported herein were 
calculated using radial increments AR equal to 0.01. 
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EVALUATION OF LOSS COEFFICIENTS ACROSS THE VANELESS FIELD 


Standard definitions of total pressure and enthalpy loss 
coefficients, in terms of the average pressures at inlet to and 
exit from the vaneless region are adopted. The total pressure 
loss coefficient is given by: 


Y 



(49) 


while the enthalpy loss coefficient may be written as: 

- 1 


P 7 ^ - 
(k^O y 


t3 


p xzi 
(_t2) Y 


(50) 


- 1 


It is clear from the above two equations that the flow conditions 
at station 3, the vane less space exit, have to be known in order to 
calculate the loss coefficients. 

In the previous chapter, it is explained how to determine, 
numerically, the Mach number, the flow angle and the effective 
nozzle depth variations along the flow path. Knowing the values 
of these variables at the nozzle exit (.station 3) , the static and 
total pressures, p^ and P^/ are computed using the equations 
derived below. 

The conservation of mass between the vaneless diffuser inlet 
and exit, is expressed as: 


p 2 V 2 COSa 2 r 2 b 2 p 3 V 3 COSa 3 r 3 I:) e 3 

where 2b is the effective depth of the nozzle at station. 3. 
e 3 

If the definition of Mach numbdr and the equation of state 

> 

are substituted into the above Equation (51) , the following 
expression for P 3 /P 2 obtained: 
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( 52 ) 


p 3 _ b 2 r 2 COSa 2 M 2 |/ 1 + y 2 ^2 
p 2 b e 3 r 3 COSa 3 M 3 V 1 + M| 


The total pressure at the vaneless space exit, P^ 3 , is then 
determined using the following relation: 



£3 

p 2 


1 + 


1 + 




Y 

Y-l 


(53) 


Equations (52) and (53) are sufficient to determine the loss 
coefficients Y and provided that the inlet flow conditions to 
the vaneless space are specified. 
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RESULTS AND DISCUSSION 


The numerical examples worked out using the present method of 
analysis are presented in two groups. In the first group, the 
effects of inlet flow conditions on the flow behavior and losses 
in the vaneless field are investigated. The second set of 
results illustrate the influence of changing the end wall spacing 
on the same flow parameters. 

1 . Effects of Inlet Flow Conditions 

The effect of two inlet flow parameters, namely the Mach number 
and the flow angle on the vaneless nozzle performance are considered 
first. Computations were carried out with an inlet Mach number of 
0.8 and different inlet flow angles ct 2 , the corresponding results 
are shown in Figures 20 to 22. On the other hand, the data obtained 
with different Mach numbers for a constant inlet flow angle of 70° 
are presented in Figures 23 to 25. The geometrical end wall spacing, 
2b, and the inlet radius to the vaneless field, r 2 , for this first 
set of results were chosen to be 0.339 and 3.93 inches respectively. 
The total temperature was held constant throughout the f low_ at 2000°R. 

The change in angular momentum ratio, r^^sina^/^v^i 11 ^ ' 
with the normalized radius, R, is shown in Figure -20 for different 
values of inlet flow angles, a 2 . It can be seen that 
increasing a 2 causes a reduction in the angular momentum ratio 
everywhere. This effect becomes particularly predominant at 
high values of inlet flow angles. For example, at a normalized 
radius of 0.9, the angular momentum ratio changes from 0.992 to 
0.987 corresponding to an increase in a 2 from 60° to 70°. 

Meanwhile, a further increase in a 2 from 70° to 80° results in 

a reduction of the same parameter from 0.987 to 0.973. 

2 2 

The radial variations of the ratio, M^/M 2 , are shown in 
Figure 2'1 for three different values of inlet flow angle. It can 
be seen from the figure that increasing a 2 reduces M 2 /M 2 at each 
radius. This is a consequence of the larger reduction in the 
tangential velocity component, due to viscosity, an influence 
that dominates the augmentation in the smaller radial velocity 

component. It can also be concluded that for all the values of a 0 

2 2 ^ 
the ratio M 3 /M 2 increases as the flow proceeds inwards which, is 

anticipated. 
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The variation in the pressure loss coefficient, Y, and 
enthalpy loss coefficient, 5 , are shown in Figures 22a and 22b 
respectively. It is clear that both coefficients increase with 
increased inlet flow angle, a 2 , at every radial location. If 
the results of Figure 22 are compared with those obtained in 
Part 1 it can be concluded that, even for the largest values of 
a 2 , the loss coefficients in the vaneless field are extremely 
low as compared to the losses in the vaned region. For example, 
the pressure loss coefficient of Figure 22b is as low as 0.0125 
at R = 0.90 for values of « 2 and M 2 of 60° and 0.8, respectively. 
On the other hand, the pressure loss coefficient of Figure 8 
is as high as 0.107 under the same operating conditions. 

Briefly, one can conclude that increasing the inlet flow 
angle, a 2 , results in a larger reduction in the total pressure, 

Mach number and angular momentum at nozzle exit. This can be 
attributed to the longer path that the flow has to take to reach 
a certain radius if the inlet flow angle is increased. 

The effect of changing the inlet Mach number, M 2 , on the 

angular momentum ratio is illustrated in Figure 23. It is 

evident that at any radius, the drop in angular momentum due to 

viscous losses, is smaller for high values of M 2 . Accordingly, 

high Mach numbers at inlet to the turbine vaneless field are 

2 2 

recommended. The variation in the ratio M 2 /M 2 is shown in 
Figure 24 for different values of inlet Mach number. It is 
obvious that the acceleration rate in vaneless field is higher 
for larger inlet Mach number. 

The pressure and enthalpy loss coefficients Y and t, are 
shown in Figures 25a and 25b respectively. It is clear that 
the pressure loss penalty paid due to viscous effects at any 
radius ratio, is not strongly affected by the inlet Mach number. 
Furthermore, the enthalpy loss coefficient at any radius remains 
practically constant for all the inlet Mach numbers investigated. 

2 . Effects of Nozzle and End Wall Spacing 

The object of the second group of numerical examples is to 
study the effect of the geometrical passage depth, 2b, on some flow 
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properties and on the loss characteristics in the vaneless 
field. For this purpose, the following values of the passage 
depth were considered: 0.262", 0.393" and 0.785" with an 

inlet radius to the vaneless field of 3.93". The corresponding 
values of the nondimensional parameter (r 2 /2b) thus ranges 
between 15 and 5. The results were obtained at 0.8 inlet 
Mach number and 70° inlet flow angle with an inlet total 
temperature of 2000°R. 

Figure 26 represents the radial variations of angular 

momentum ratio, r^V^sina-j/^V^inc^ , for the three values 

of the parameter (r 2 /2b) . Although the angular momentum 

ratio decreases at any radius, R,'with the decrease in end 

wall spacing, 2b, the curves indicate that there is only minor 

differences between smallest and largest values of the 

2 2 

parameter (r 9 /2b) . The change in the ratio, M 2 /M 2 , with R 

is shown in Figure 27. Increasing the passage depth, 2b, 

2 2 

results in an increase in M 9 /M 9 at any radius ratio, R. 

2 2 J z 

The increase in M 2 /M 2 results primarily from the reductxon 
in the friction surface area compared to the total flow area 
as the wall spacing is increased. Figure 27 also indicates 
that the variation of nozzle end wall spacing, 2b, has a 
small influence on the Mach number distribution within the 
vaneless space. 

The influence of the end wall spacing on flow losses is 
shown in Figures 28a and 28b. It is clear that larger losses 
are incurred as the nozzle passage becomes narrower, corresponding 
to high values of r 2 /2b. It is important to emphasize once 
more that, even for the large values of (r 2 /2b) , the losses 
in the vaneless field remain relatively low as compared to 
the vaned region losses. 
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CONCLUSIONS 


An analytical model which is based on the observations 
obtained from an experimental study was developed to predict the 
losses in a radial nozzle annulus. The contributions of the 
end wall boundary layers to the losses was found to be several 
times that of the vane surface boundary layers. This influence 
is particularly significant when the end wall boundary layers 
are characterized by large cross flow components. Under the 
stipulation of the analytical model in which the mixing losses 
are lumped with the friction losses in the vaned region, a small 
portion of the overall losses were contributed by the viscous 
effect in the vaneless field. Experimental findings in which 
stator-rotor interaction effects are not considered show that 
the losses are influenced by flow nonuniformities at inlet to 
the nozzle channel and by secondary flow. These factors need 
further investigation if a realistic estimate for the losses 
is to be obtained analytically. 

It is concluded that, generally the losses in a radial 
nozzle assembly would not be greatly affected by the addition of 
a large vaneless space. Also, if the nonuniformities of the 
flow discharged from the vaned region can be reduced, the 
efficiency of the assembly as a whole would be improved. Thus, 
the loss penalty paid during an expansion process, resulting in 
the required flow properties at rotor tip from a specified 
outlet flow condition from the scroll, could be minimized. 

Such minimization is achieved by a proper selection of vanes 
configurations in conjunction with a suitable vaneless space 
dimension. The application of the compositive method of analysis 
presented allows one to differentiate between various suggested 
radial nozzle designs in order to select the optimum configuration. ' 
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NOMENCLATURE 


Symbol 

A 

B 

2b 

2b e 

C f 

G 

g 


K, L 


M 

M* 


n 

P 


P 


P 


t 


r 


R 



T 


V 


parameter (see Eq. 25b) 
depth ratio defined as 2b & /r 2 

axial depth of the no 2 zle - distance between the 
two end walls , ft. 

effective depth of the nozzle 

skin friction coefficient 

streamwise profile function 

cross-flow profile function (see Eq. A7) 

streamwise shape factor for the end walls boundary 
layers (see Eq. A8) 

profile parameter, for the end walls boundary 
layer (see Eq. A9, A10) 

Mach number, dimensionless 

critical Mach number, dimensionless 

coordinate normal to the vane surface, ft. 

pressure ratio defined as P/P t2 

2 

static pressure, lb f /ft 

2 

total pressure, lb^/ft 
radial coordinate, ft. 
radius ratio defined as r/r 2 

Reynolds number based on momentum thickness of 
the main flow, dimensionless 

temperature, °R 

vane trailing edge thickness, ft. 

streamwise velocity component in the boundary 
layer over the end walls , ft/sec 
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V 


velocity, ft/sec 

w cross-flow velocity component in the boundary 

layer -over the end walls, ft/sec 

Y overall total pressure loss coefficient 

z axial coordinate, ft. 

Z number of vanes, dimensionless 

a angle between flow direction and inward radial 

direction, radians 

8 main flow turning angle, radians (see Fig. 14) 

y ratio of specific heats, dimensionless 

T parameter (see Eq. 6A) 

5 boundary layer displacement thickness, ft. 

trailing edge blocking factor defined as 
te/— ^ r^cosa-^) 

* 

6 total nondimens ional displacement thxckness over the 

2 ^ 

vane surfaces defined as (6 s +6p)/{— r^cosa^) 

"it 

6 streamwise total nondimens ional displacement thickness, 

over the two end walls, defined as (<5 UX + 6^ x )/2b 

9 boundary layer momentum thickness, ft.; or angular 

' coordinate in plane normal to axis of rotation, radians 

* 

0 total nondimens ional momentum thickness over the vane 

2 TT 

surfaces, defined as (6^+6 )/(— g- r^cosa-^) 

6 streamwise total nondimens ional momentum thickness 

over the two end walls, defined as C 0 UX +0 £ X )/2b 

X parameter (see Eq. A23) 

p density, lb/ft 2 3 

n parameter (see Eq. A31) 

ij> w angle of limiting wall streamline and with 

respect to V fs 
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x Q shear stress at the wall 

5 overall enthalpy loss coefficient 

s tangent of the limiting streamline defined as tanij) 

Subscripts 

0 inlet to radial nozzle, outlet from the scroll 

1 nozzle exit 

/ 

2 station at which uniform conditions are assumed to 

/ 

take place 

3 inlet to rotor 

c cross flow direction 

cr critical conditions 

« local conditions outside the boundary layer regions 

6 end wall lower surface 

p vane pressure surface 

s vane suction surface 

t total conditions 

u end wall upper surface 

r radial direction 

x streamwise direction 
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APPENDIX A 


FORMULATION OF RELATIONS BETWEEN INTEGRAL QUANTLTIES 
AND BOUNDARY LAYER CHARACTERISTIC PARAMETERS 


A.l Basic Relations and Definitions 


In this Appendix, the integrals which appear in the left 
hand side of Egs. (2) , (4) and {6) will be determined in terms 
of the boundary layer characteristics at station 1. To 
accomplish this, the domain of integral will be divided into 
a main stream region and boundary region as shown in Figure Al. 
This is a generalization of the approach used in the two 
dimensional analysis of Reference [5] . 

At exit from the nozzle channels the profile boundary 
layer developed along the vane surfaces is collateral, while 
the end wall boundary layer is skewed. The velocity vectors 
within the skewed boundary layer can hence be resolved into a 
streamwise component v, and a cross flow component w. The 
displacement and momentum thicknesses are used to describe the 
boundary layers at station 1. These two quantities for the 
collateral boundary layer are defined as : 


5 


5 * - 


/ 

o 


PV 


TpvT 


dn 


(Al) 


e 



pV 


(PV) 


dn 


PV" 


(pV 2 ) 


dn 


(A2) 


Where n is the coordinate normal to the vane surface and <5^ 
is the boundary layer thickness. On the other hand, the streamwise 
and cross flow end wall boundary layer displacement and the 
momentum thickness are defined as follows: 


x 


= 6 ^ - 


f 

r pv 

i w 


dz 


(A3) 
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x 




S 


2 

f pv 


(pV 2 ) 


dz 


CO 


5 = 


f 

r p w 

L w 


dz 


B - f 


pw 


dz 


o (pV‘) 


(A4) 


(A5) 


(A6) 


Where z is the coordinate normal to the end wall surface. 

The relation between the parameters of equations (A3) through 
(A6) can be determined for a given boundary layer profile. 

It is assumed that in the end wall boundary layer the streamwise 
and the cross flow velocity profiles are given by Reference' [6] . 



G 



eG g 


(A7) 


In the above equations G and g are functions of the coordinate 
perpendicular to the end wall/ z, normalized with respect to the 
boundary layer thickness (<5f) , "s is the tangent of the angle, 
between the surface shear and the local main flow direction 
as shown in Fig. A2. 

The following streamwise and cross flow parameters are 
used to describe the end wall boundary layers using the assumed 
formes of the velocity profiles: 


H x = e 


X 


X 


r 

o_ 

f 

o 


(1-pG) dz 


p{l-G)G dz 


(A8) 
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K 



L 


0 

c 


e 2 0x 


/ y(Gg)dz 
o 

6 f 

/ (y-yG)G dz 
o 


/ y (Gg) dz 
o 


f (y-yG) (G)dz 
o 


(A9) 


(A10) 


Where y is the nondiraensional ratio p/p m . 

In Ref. [7] various families of functions G and g that are 
most compatible with the experimental data were determined. 

The corresponding profiles are plotted in Fig. A3, together 
with the Band, "B" , of the experimental data obtained by 
Johnson [6]. The various parameters of Eqs. (A12) , (A13) and 
(A14) were calculated for the incompressible flow corresponding 


to the profiles of 
table . 

Fig. A3 

, and are given 

in the following 

Profile 

H 

X 

K 

L 

A 

1.286 

0.457 

0.0359 

B 

1.37 

2.43 

0.968 

C 

1.40 

2.27 

0.989 

D 

1.286 

1.249 

0.262 

TABLE A-l : TYPICAL 

VALUES 

FOR THE END WALL BOUNDARY LAYER 


CHARACTERISTIC PARAMETERS. 
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A. 2 Evaluation of the Different Integral Terms in the 
Governing Equations 

The integrals in the Governing Equations (2) , (4) and (6) 
will be expressed in terms of the boundary layer parameters of 
Eqs. (Al) through (AlO) based on the following assumptions: 

1. Different boundary layers do not interfere with 
each other, hence corner effects are neglected. 

2. The flow is homogeneous outside the boundary 
layer regions at exit from the nozzle channels. 

In accordance with the foregoing assumption, the velocity 
distribution of Figure A3 due to the profile boundary layers 
on the vane surfaces will result in the variation of the total 
pressure in the tangential direction shown in Figure A4. The 
total pressure and velocity deficiencies in, the axial direction 
caused by the end wall boundary layers will be similar to those 
indicated in Figures A5 and A6 . 

Referring to Figure A2 , the following equations relate 
the different skewed end wall boundary layer velocity components 
and flow angles: 


(A12) 
(A13) 

and V 2 cos 2 a may be approximated by 


V cosa = v cosa^ + w sina-^ 

2 2 . 2 . 

V cosa sina - v cosa^ srna^ - w cosa^ sma- 


2 2 
V cos a 


2 2 ^ 2.2 
v cos + w sm 


(A14) 


Where a is the angle between the flow direction and the radial 
inward direction as before, while < 3 ^ is the same angle for 
the local free stream at station 1, which is independent of z 
and 8 according to the assumptions. 
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+b -Ht/Z 

A. 2.1 Evaluation of the Integral f [ / (p Vr cosa) d 0 ] dz 

-b -ir/Z 

Referring to Figure Al, it is more convenient to breakup 
the above integration into three separate integral's as follows: 

+b +ir/Z +b-<S u f +*/Z 

/ [ / (pVr cosa)..de]dz = cosa.r., f If {pV),d 8 ]dz 

-b -ir/Z x X -b+5^ _7r / z 1 

" b+<S & f +u/Z b +ir/Z 

+ r-, { / [ / (pV cosa) . d0 ] dz + / [ / (pV cosa),d9]dz} 

-b -tt/Z 1 b-S.. -ir/Z 

£ (A15) 

where 

6 , 6 are the full boundary layer thicknesses over 

Uf * £ 

the upper and lower end wall surfaces. 

Referring to Figure A2 , the coordinate used in the profile 
boundary layer equations is related to the circumferential 
coordinate 0 by: 


de 


dn 

r^cosa 1 


(A16) 


Substituting the above expression for d 8 in the first term in 
the right hand side of equation (A15) , and using the definition 
of the displacement thickness of Equation (Al) we can arrive 
at: 


+b- 6 1 

, U f +ir/Z 

J [ / (pV),de]dz = [2b -(«„ + 6 )]{ pV) 

-b+S 0 -tt/Z X U f *£ 

*f 


'TT 


Z ^ {1 " 6 5 te ] 


(A17) 


where 


5 * - (5 S + y/f 1 r i cosll i 


(A18) 


2-jt 

5 te = te/ Z~ r l COSct l 


(A19) 
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and 


d s , 6 are the boundary layer displacement thickness 

over suction and pressure surfaces of the vane. 

te is the vane trailing edge thickness. 

2 77 , t 

— — . is the angular spacing between two successive 

vanes. 

Using Eq. (A12 ) , the last two terms on the right hand side of 
Equation (A15) could be rewritten as: 


-b+<S 

/ 

-b 


'f +ir/Z +b 4-ir/Z 

[ / (pV cosaKdeJdz + / [ f (pV cosa)-,d9]dz 

-tt/Z +b-<S ~tt/Z 




-b+6 a ^ p (v + w tana.) 

z r cpv) » 1 cos «i t / — rtvjz d2 


+b 

+ / 
+b-6 


p <v + w tana^) 

(pvTT 


dz] 




When Equations (A3) and (A5) are used, the above equation 
reduces to : 


/ 

-b 


f 

+ tt/Z 

+b 


E f (pV coset) . de] dz 

+ / 


-tt/Z 

+b-6 

2tt 

z 

(pV)^ r l COSct l E + 

5 u 


+ir/Z 

t / (pV coses) .dejdz 

/r-7 -L 




where 


5 £x' 5 ux 


<S„ . 5 

£c' uc 


are the streamwise boundary layer displacement 
thicknesses over the lower and upper end walls 
respectively. 

are the cross flow boundary layer displacement 
thicknesses over the lower and upper end walls 
respectively. 
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Finally substituting Equations (A17) and (A20) into (A15) 
we get: 


r2ir, 


+b tt/Z 

/ [ / (pVr cosa^deldz = (pV) M r ] .cosa 1 [-2b- {-6 U +6^ )] 

-b -n/Z f f 


*[1-6 - fi te 3 + [5 Jl f +,S u f ~ (6 £x +5 ux )+C5 £c +6 uc )tana l 3 


(A21) 


If we define the following non dimensional parameters 


** 


0 ux + 0 £x 
2b 


(A22) 


A - H - eK tana. 


(A23) 


where 

9 r 0 are the streamwise boundary layer momentum 

ULK 36 X 

thicknesses over the upper and lower end 
walls respectively. 


Equation (A21) can be simplified by neglecting the higher 

* 

order terms (6 +6 )6 , (6 -+-6 ) 6 . and introducing 

u f z f u f z f te 

the parameters defined by equations (A8 ) , (A9) , (A22) and (A23) 

giving the following relation: 


+b +ir /Z 

/ [ / (pVr cosa),de]dz = (pV) 

-b -ir/Z ~ i 


2 IT 
Z 


r^cosa^ 


2b[l-6*-6 te -X6**] 

(A24) 


+b +tt/Z „ 2 

A. 2. 2 Evaluation of the Integral f [ f (r pV cosasina) , d8 ]dz 

-b —tt/Z 

The evaluation of the integral will be divided over the 
three same regions discussed previously. 
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+b +ir/Z 2 2 2 

/ [ / (r pV cosasina) ,de]dz = {r, cosa^sina,} [2b - 

-b -ir/Z X 111 

+ir/Z 2 2u 2 “ b+6jl f 2 

-(<$, +5 )] / (pV ),de + (■=— r,) [ / (pV cosasina) . dz 

U f — ir/Z 1 Z 1 ~b 1 


b 

+ / 


b-6 


(pV cosasina) ^dz] 


(A25) 


u, 


Using Eqs. (Al) , (A2) and (A16) into the first term in 

the right hand side of Eg. (A25) , it reduces to: 

2 +*/Z 2 

{r, cosa, sina, } [2b - (5 + 6 ) ] f (pV~),de 

U f *f -ir/Z 1 


“ (PV 2 )^ sina 1 cosa 1 t2b-(6 u +6 £ )]|^rj [1-6 *-6^-0*] (A26) 


where 

0 = (8 g + 9 ) /g— ’ r 1 cosa^ (A27) 

and 

9 and 9 are the boundary layer momentum thickness on 

° ir 

the vane suction and pressure surfaces 
respectively. 

When Eqs. (A13) , (A4) and (A6) are substituted into the second 

and third tern of (A25) , it can be deduced that: 

-b+6 £ ^ 

r 2 2 

f , {cosasina pV ) ± dz = (pV ) w cos^sina^fi f -8 A -e ax **0 3 

— b 1 

(A28) 
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and 


b 2 
/ (cosasina pV ),dz 

b-6 

u^ 


(pV 2 ) 


cosa,sina, [6 _-<$ -0 -6 j 

1 1 L uf ux ux uc J 


(A29) 


where 


, 8 and 0 are the cross flow boundary layer momentum 

uc Jo c 

thickness over the upper and lower end walls , 
respectively . 

Finally, substituting Equations (A26) , (A28) and (A29) into 

(A25) and making use of the parameters given by (AS) and (A10) 
the integral of equation (4) is expressed as follows: 

+b +ir/Z 2 2 2 2ir 2 

/ [ / (r pV cosasina) ^d8] dz = (pV ) ^ sina^cosa^ ^ — r^ 2b 

-b -tt/Z 1 

£ Jc icit 

* [ 1-6 - « te - 0 - ne ] (A30) 


where 

n = 1 + H + e 2 L (A31) 

A 


+b +ir/Z 2 

A. 2. 3 Evaluation of the Integral I [ S (pV cosa rcos (a+0) 1 d8] dz 

-b -tt/Z 

Using trigometric relations the above integral may be 
expanded as follows : 

+b +ir/Z 2 +b +ir / Z 2 2 

f [ f (pV r cosacos (a+8) , de ] dz = /[ / (pV r cos a) , cos 8d0 ] dz 

-b -ir/Z 1 -b -tt/Z 

+b +ir/Z 2 

- / [ / (p'v r cosasina), sin8d0]dz (A32) 

-b -ir/Z X 

The first integral in the right hand side of Equation (A3 2) 
may be rewritten as : 
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+fo +ir/Z 2 2 2 

f [ / (pV r cos a) cose d0]dz = {r,cos o,} -U2b - <$ 

-b — ir/Z 1 1 1 


+tt/Z +T7/Z *' b+6 & f . ? 

-5 ] / (pV ) 1 cos8d8} + r, { / [ / (pV cos a) , cosedz] de 

f -tt/Z -tr/Z -b 1 


+ir/Z b „ 

+ / [/ (pv cosa ) _ cos9dz]dS} {A33) 

-tt/Z b-5 1 


For unseparated flows, the angular coordinate, 9, does 
not change appreciably inside the boundary layers formed over 
the vane surfaces. Consequently, 9, can be considered equal 
to the angular spacing, ir/Z, within these boundary layer 
regions. Using this value of 9 together with Equation (A2) , 
it can be easily deduced that: 

+tt/Z ? ? 9 + 0 +ir/Z 

/ (pv Kcosede = -(pV ) cos ~ <-■■■■■ p ) + v- f (pV),cos9d9 
-i/Z 1 "1 2 r l cosa l fs l -,/2 1 

(A3 4) 


Using the definition of the displacement thickness given by Equation 
(Al) , the second term in the right hand side of Equation (A34) 
reduces to: 


+tt/Z 

f 

-tt/Z 


(pV) 1 cos0d9 = 



r sin tt/Z 
1 tt/Z 


COS |r (6 


5 te>' 


(A35) 


If Equation (A35) is substituted into (A34) the following relation 
can be written: 


+ir/Z 2 

/ (pV ) cos9d9 = 
-tt/Z x 



r sin tt/Z 
1 tt/Z 


cos y + 5 te + 0 ^ 


(A3 6) 


Using Equations (A14) , (A4) and (A6) , the integrals in the right 
hand side of (A3 3) reduce to: 
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+ir/Z 

f 

-ir/Z 


-b+5 , 

/ 

~b 


'f 2 2 2 2 

(pV cos a) 1 cos0 dzde = 2sin £ cos a, (pVl (5 
x z x do ^ j 


- 5 <i - e n + \ tan “1> 

X X c 


f 

(A37) 


tir/Z +b 
f f 
— ir/Z +b-<5 


(pV cos aKcosedzde - 2sin \ cos 2 a, (pV 2 ) (5. 

J- Z X 03^ 


U, 


■ 6 U - 9 U + 0 U tan “l> 
X X u c 


(A38) 


Substituting Equations (A36) , (A37) and (A38) into the right hand 
side of (A 33} , we arrive at: 


+b tir/Z 2 2 

/ [ / (pV r cos a) cosedQ 3 dz = 2b 

-b -tt/Z X 


2tt 

Z 


r^cos 


(pV 2 ), 


. 7T 

sxn z 

jr_ 

Z 


[1 


** 2 2 -it * 

- 0 (1 + - s L tan ct-j^) 3 “ cos (6 + 


te 


* 

+ 0 ) 


(A39) 


2 

where H , e L are given by the Equations (A8 ) and (A10) . 
Following a procedure similar to the one used to obtain equation 
(A39) / it may be easily deduced that: 

+b +ir/Z 2 

/ [ / (pV r cosasina) , sinede} dz = 0 (A^O) 

-b -u/Z 1 

Finally, substituting Equations (A39) , (A40) into (A33) , the 

required value of the integral in Equation (6) is obtained as: 


+b +ir/Z 2 

/ / [pV cosa rcos (a+0) 3 n dedz = 2b(pV 2 ) cosa., ~ r, 

-b -ir/Z 1 “1 1 z 1 


TT 

s in "f? ^ 

1-jZ U -re* )- 

Z 


cos Ju + 6 te + 0*) ] 


(A41) 
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where 


r = 1 + - e 2 L tan 2 a 1 (A42) 

Equations (A24) , (A30) and (A41) will be used in equations 
(2) , (4) and (6) to express the continuity, the angular momentum 
and the linear momentum equations in terms of the boundary 
layer characteristic parameters at station 1. 
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FIG A-l FLOW REGIONS AT EXIT FROM THE NOZZLE CHANNEL. 



MAIN FLOW STREAMLINE 



FIG, A- 2 NOMENCLATURE AND SYMBOLS AT VA N E S - TRAI L I N G EDGE 
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fig:a-4: typical velocity distribution in 

CIRCUMFERETIAL DIRECTION AT 
NOZZLE EXIT 



FIG. A- 5: TYPICAL PRESSURE DISTRIBUTION IN 

CIRCUMFERENTIAL DIRECTION AT 
NOZZLE EXIT, 
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UPPER SURFACE 


LOWER SURFACE 



FIG.A-6: TYPICAL TOTAL PRESSURE Df STRIBUTION 


IN AXIAL DIRECTION AT NOZZLE EXIT. 
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FIG. A- 7 TYPICAL VELOCITY DISTRIBUTION 


AXIAL DIRECTION AT NOZZLE EXIT. 
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